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ABSTRACT 

Observations from the Heliospheric Imagers (HI-1) on both the STEREO spacecraft have 
been analysed to search for very long period large amplitude stellar variability, finding 6 new 
candidates. A total of 85 objects, mostly previously known Mira variables, were found to 
show convincing variability on time scales of over a hundred days. These objects range in 
peak brightness from about 4th magnitude to 10th magnitude in R and have periods between 
about 170 days and 490 days. There is a period gap between 200 and 300 days where no 
objects were found and this is discussed. 15 of the Miras in the sample are previously recorded 
as having variable periods and the possibility for these and 2 other stars to have undergone 
a period change or to be irregular is discussed. In addition to the 6 stars in the sample not 
previously recorded as variable, another 7 are recorded as variable but with no classification. 
Our period determination is the first to be made for 19 of these 85 stars. The sample represents 
a set of very long period variables that would be challenging to monitor from the Earth, or 
even from Earth orbit, owing to their position on the Ecliptic Plane and that their periods are 
often close to a year or an integer fraction thereof. The possibility for the new candidates to 
possess circumstellar shells is discussed. 
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1 INTRODUCTION 

Roughly four hundred years ago the first variable star was recog- 
nised due to its changing brightness (Hoffleit 19971. This red ob- 
ject, showing striking light variations in the visual, was later named 
o Cet and became the first member of the variability class called 
long period variables (LPVs) today. Being easily detectable be- 
cause of the large photometric amplitudes they were intensively 
studied thereafter and represent now a prominent group within 
the General Catalogue of Variable Stars (Samus et al. 2012). It is 
known nowadays that LPVs are stars of low to intermediate main 
sequence mass (~0.8-8M fJ ) in a quite late stage of stellar evo- 
lution, the Asymptotic Giant Branch (AGB). Such late-type gi- 
ants populate regions of high luminosities (several Lq) and 
low effective temperatures (below ~4000K) in the Hertzsprung- 
Russell diagram. During the AGB evolution, the stars start to pul- 
sate which causes the pronounced photometric variations (e.g., 
|Qlivier & W ood 20051. Conventionally, LPVs were subclassified 
into Miras, semiregular variables, and irregular variables with the 
visual amplitude and the regularity of the lightcurve as criteria. An 
important step forward in our understanding of LPVs was provided 



by the large surveys of the Magellanic Clouds and the Galactic 
Bulge (e.g., |Wood||2000| |Ita et al.||2004| |Groenewegen & Blom-| 



|maert||2005| |Matsunaga, Fukushi & Nakada| 2005|>, recently up- 
dated with space-based IR data i Riebel et al.||2010| >. From these 
we know that AGB variables constitute a few period-luminosity re- 
lations. The most prominent LPVs, namely the Mira variables, can 
be found along the sequence for fundamental mode pulsators. The 
majority show large-amplitude variations (AV of a few magnitudes, 
AK £3 0.4-l ma S) on time-scales of 100-600 days (e.g., |White-| 
|lock, Marang & F east 2000). The P-L relation for Miras was de- 
termined rather precisely based on a number of well-characterised 
LMC objects by Whitelock, Feast & van Leeuwen ( 2008} but could 
be investigated also in other Local Group Galaxies (e.g., Lorenz et| 
|al.|2011) and even in systems as distan t as Cen A (|Rejkuba|2004[ >. 
Apart from pronounced mass loss (e.g., [Nowotny et al. [201 \) AGB 
stars are also characterised by the occurence of He-shell flashes 
(thermal pulses; Lattanzio & Wood 2004). Not only are these re- 
sponsible for changing the atmospheric chemistry from O-rich to 
C-rich via dredge-up of carbon from the interior, thermal pulses 
are also suspected to be the reason for period changes in observed 
lightcurves (e.g., Templeton, Mattei & Willson|2005| |Uttenthaler| 



|et al.|201 H|Lebzelter & Andronche|201 1\ . 



* e-mail: k.t.wraight@open.ac.uk 



The data used in this study comes from NASA's STEREO mis- 
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sion, which aims to image the Sun's corona in 3D and observe 
Coronal Mass Ejections (CMEs) from the surface of the Sun out 
to the Earth's orbit ( [Kaiser et al.||2008} . The two satellites are in 
different heliocentric orbits, one slightly inside the Earth's orbit 
(STEREO-Ahead) and one slightly outside (STEREO-Behind) and 
the angle between each satellite, the Sun and the Earth increases 
by about 22.5 degrees every year. Photometry of background stars 
in the images is possible as the calibration of the Heliospheric Im- 
ager cameras has been performed to a very hi gh standard jBrown 
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Bewsher & Eyles|2009[ |Bewsher et al.|2010| [Bewsher, Brown& 



Eyles|2012||Wraight et al.|201 1| 2012). Having two nearly identi- 
cal satellites in different heliocentric orbits provides greater phase 
coverage of long period variables, with greater homogeneity than 
achievable from observations conducted from the ground. 

In the following sections we outline the characteristics of the 
STEREOfHl-l observations and explain how the sample was ex- 
tracted and analysed. We provide periods for 85 stars showing 
large amplitude variability 0.3 magnitudes) with periods longer 
than 100 days. The majority of these are known to be Miras or 
semi-regular variables. We estimate times of maximum brightness, 
where visible. Our analysis of patterns found in the sample focuses 
on the errors obtained for the periods and also on a gap in the pe- 
riods found between 200 and 300 days. We discuss a number of 
individual stars, in particular those which are known to have vary- 
ing periods, those which have not previously been classified and 
those found here to be variable for the first time. 



2 OBSERVATIONS AND ANALYSIS 

2.1 Characteristics of STEREOIWl-l observations 



The Heliospheric Imagers are described in detail in Eyle s et aL| 
(2009 ) but we summarise the main features here for convenience. 
The field of view of the STEREO/Hl-l cameras is 20 degrees by 20 
degrees, centred 14 degrees away from the Sun's centre. The aper- 
ture of the cameras is just 16 mm and the focal length is 78 mm. 
Each CCD is 2048x2048 pixels, which provides a plate scale of 35 
arc-seconds per pixel. The images are binned 2x2 on-board, to re- 
duce the bandwidth, and as a result the final images received have 
a resolution of 70 arc-seconds per pixel. Each image is the result of 
30 summed exposures, with each exposure lasting 40 seconds, and 
one image is produced every 40 minutes. In this way, stars as faint 
as 12th magnitude in R can be observed but only the very bright- 
est stars will saturate the CCDs. Stars of about 4th magnitude may 
show some systematic effects due to saturation, whilst 3rd magni- 
tude stars often show systematic effects but may be usable whereas 
stars of 2nd magnitude and brighter are unusable. We use the NO- 
MAD 1 catalogue (Zacharias et al. 2004 ) to determine which stars 
to observe, selecting all those listed as 12th magnitude or brighter 
in the R band within our field of view. Aperture photometry is per- 
formed on each star as described in Bewsh er et al.| ( |20 10) and the 
data reduction pipeline is summarised in Wraight et al. (201 lj. This 
provides an extremely useful resource, with photometry of almost 
900,000 stars along the Ecliptic Plane. 

The STEREO/HI- 1 imagers have an unusual spectral bandpass 
(Figure [TJ. They are most sensitive inbetween 630 nm and 730 nm 
but also have some sensitivity in the blue, around 400 nm, and in 
the infra-red at around 950 nm. The sensitivity in the blue is very 
useful for observing hot stars (Wraight et al. 2012) but it is the sen- 
sitivity in the infra-red that makes the observations of many of the 
very cool stars in this paper possible. Some of the stars in our sam- 
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Figure 1. Plot of the filter throughput convolved with the quantum 
efficiency of the CCDs for STEREO/Hl-iA, shown in a dashed line, 
(STEREOfHl-lB is almost identical (Bewsher et al. 2010), shown here in 
a dot-dashed line). This is set against a synthetic spectral energy distribu- 
tion of a red giant with MK spectral type M7 (Fluks et al. 1994), shown in 
a dotted line, with the continuum emission in a solid line. 



pie are cool enough, or sufficiently obscured by circumstellar ma- 
terial, to be fainter than 12th magnitude in the R band and do not 
feature in the STEREO database directly. In these cases the large 
pixel size of the STEREOfHl-l imagers, combined with the sensi- 
tivity near 950 nm, allows them to be observed indirectly through 
blending with a nearby star that is in the database. As Mira vari- 
ables are known to change their temperature and spectral type, the 
magnitude of the variability observed by STEREOfHI-l is not reli- 
able. This can be seen by referring to synthetic spectra of red giants 
( |Fluks et al.|[T994> . There is a significant difference in the emis- 
sion in the region 630 nm to 730 nm between the spectral type M6, 
where there is still some emission, and M7 and later, where there is 
relatively little emission. A star varying across this threshold will 
therefore have an exaggerated amplitude. Giants of spectral type 
M7 or later are seen almost entirely through their emission in the 
infra-red and may show different variability to that seen in other 
parts of the spectrum, whilst giants of M6 and earlier may show 
variability due to processes affecting both the red and the infra-red. 
In particular, it is worth noting that some of the faintest stars in 
our sample, including many of the new candidates, are so faint that 
only deep all-sky surveys e.g. Zacharias et al. ( 2004 ) and Monet 
|et al.| p003 l have detected them in visible light. Even in the infra- 
red observations are scarce, mostly fro m the IRAS mission jNeuge-| 
|bauer et al.|1984| >, the AKARI mission (Murakar ni al.|2007) and in 
the 2-Micron All-Sky Survey (2MASS) ( |Skrutskie et al.|2006| >. 
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2.2 Quality of STEREO/H1-1 data 

In order to determine a period for the candidate variables, it was 
necessary to use undetrended data and the only processing that 
was applied was the exclusion of obvious outliers and system- 
atics, where possible. This is because the polynomial detrending 
normally used to clean the lightcurves removes long period vari- 
ability. As a result, many of the lightcurves show indications of 
artificial trends, often a result of the flat-fielding breaking down 
near the edges of the detectors. Other sources of noise affect the 
long period signals being searched for to a lesser degree, most no- 
tably de-pointing events associated with micrometeorite hits, which 
are more common in the data from the STEREO-Behind satellite, 
STEREO/Hl-IB, also observed by |Davis etal] p012) >. The pres- 
ence of artificial trends and the ability to recognise them as such 
also limits the number of maxima that can be reliably observed for 
the stars in the sample and might also influence the times of max- 
imum light, in those cases where the effect is small enough to be 
reasonably sure a maximum is genuinely being observed. 

2.3 Data analysis 

In order to extract the sample, the first step was to select those stars 
showing the largest difference in the weighted mean magnitude ob- 
served by the two satellites. We anticipated that the majority of 
these would be due to systematic effects, mostly relating to flat- 
fielding near the edges of the CCDs but also due to planetary incur- 
sions from Venus and Mercury which frequently pass through the 
field of view. The photometry is also conducted very slightly dif- 
ferently between the two imagers, using aperture photometry with 
different apertures iBewshe r et al.|20l0| >. As a result, there are dif- 
ferences between the weighted mean magnitudes observed by the 
two satellites, with these systematics being greatest in the Galactic 
centre and anti-centre. In normal circumstances, however, the mag- 
nitude of these differences amounts to less than 50 mmag and ex- 
perience has found empirically that systematics are not significant 
unless the magnitude of the difference between the two satellites 
is larger than 0.1 mag. In contrast, the smallest such difference of 
any object included in the sample is that of V901 Sco, a known 
semi-regular variable, which in STEREO/Hl-l shows a difference 
of slightly over 0.3 mag between the two satellites and is unlikely to 
be due solely to systematic effects. The presence of these system- 
atics is the reason why Miras especially were searched for, as their 
variability is so large with respect to this background noise that 
mis-identifications are relatively unlikely. Approximately 10,000 
lightcurves were visually examined showing large differences in 
the weighted mean magnitudes between the two satellites, from 
which a sample of about 130 was recovered that appeared to be 
more likely due to genuine variability than any known systematic 
effects. The possibility of this process introducing a selection ef- 
fect is discussed after the results have been presented, as it is one 
possible cause for a lack of any objects in the sample showing a 
periodicity between 200 and 300 days. For every one of these stars, 
some basic information was extracted from the SlMBAD database 
and the NOMAD 1 catalogue (Zacharias et al. 2004}, so as to iden- 
tify the source of the variability observed. In many cases, espe- 
cially in the Galactic centre, it was not possible to be reasonably 
sure of the origin of the observed variability, mostly because there 
were too many candidates but often because there were no can- 
didates showing the colours expected of a red giant. A few were 
deselected after further examination showed that systematics were 
more likely responsible for the variability, or because a period could 



















- 














- 

• 


• • 












« 

• 


• 


••• 


• 

• 








-» — 






■ 






— •- 

• 










• i 


t *iC - 













































50 100 150 200 250 300 

RA (DEGREES) 



Figure 2. Plot of the locations of the 85 stars in the sample on the sky, 
with Right Ascension and Declination given in degrees. As STEREOfUl-l 
only observes stars within 10 degrees of the Ecliptic Plane, this is where 
the candidates are found. Note that more are observed near to the Galactic 
Centre. 



not be found during the final stage of the analysis. A sample of 85 
was eventually analysed in detail using PERANSO Q and a period 
determination made using a Discrete Fourier Transform ( Deeming 
1975|> (he reafter DFT), Phase Dispersion Minimisation ( Stellingw- 
erf| 1978[ > (hereafter PDM) and the Renson string length minimisa- 
tion method l |Renson|19 78 ) (hereafter SLM). Although PERANSO 
has numerous algorithms implemented, these three have a funda- 
mentally different basis to each other and are expected to have dif- 
ferent, hopefully complementary, strengths and weaknesses. It was 
important for the methods chosen to be able to deal with a very 
small number of data points, as each epoch of about 20 days would 
effectively be treated like a single point. The ability to straight- 
forwardly and manually remove known, obvious, artefacts, such as 
planetary incursions due to Venus and Mercury passing through the 
field of view, made PERANSO preferable to other programs, at the 
expense of being slightly more time-consuming, which was why 
using numerous additional algorithms was not done. The median 
value found by the three algorithms was used, so that if any two 
algorithms agreed this would produce a more reliable period. This 
produced a proportional difference between the periods found and 
the previously known periods of about 4% (Section |3 - 1 1 > - Periods 
were searched for in the range of 50 to 1000 days and the strongest 
signal of each of the three different algorithms was recorded, along 
with the error. Where possible, times of maxima were also recorded 
along with the corresponding magnitude observed. These periods 
and ephemerides are one of the main results of the research pre- 
sented herein, along with the photometry presented in the form of 
lightcurves phase-folded on these periods. The location of the 85 
stars in the sample on the sky is shown in Figure [2] As an example 
of a bright well-known Mira, we show the undetrended data and 
the phase-folded lightcurve for R Cnc along with the periodogram 
from a DFT analysis of this data in Figure[3] 



http://www.peranso.com 



4 K. T. Wraight, D. Bewsher, Glenn J. White, W. Nowotny, A. J, Norton and C. Paladini 




Figure 3. The known Mira variable R Cnc, as observed by STEREO/HI- 1 . The raw lightcurve is shown on the left, the DFT of this lightcurve after culling of 
outliers is then shown in the middle plot and the phase-folded lightcurve on the strongest period (179 days) is shown on the right. Note that the known period 
of this well-observed star is 361.6 days I Samus et al. 2012 1. 



3 RESULTS 

The analysis produced periods, with errors, for 85 stars and the 
time of at least one maximum for 24 of these (Table [TJ but V932 
Sco was excluded from the table as it is a known Orior0 variable 
and the table is likely to be used as a resource for bright LPV stars. 
This is the first period determination for 19 of these 85 stars, 7 of 
which were previously unclassified and 6 previously unknown to be 
variable. Table [2] shows when the 6 previously unknown candidate 
variables will next be observed by STEREO/HI- 1 . The phase-folded 
lightcurves of all 85 stars on the periods found are given in the 
Appendix (Figures |A~T) | A2 1 1 A3 1 1 A4 1 and | A5 } . Unfolded lightcurves 
for some stars of interest are also shown (see Section[4j. V932 Sco 
has no period given in |Samus et al.| l[2012) and the period found 
here of 423 ± 87 days is the first determination of a period, however 
no maxima were observed. As it is important to remind the reader 
that there is a chance that the new variables (marked with " — " in 
column 1 1 of Table[T]( and the unclassified variables (marked as V* 
in column 1 1 of Table[TJ are not necessarily Miras or semi-regular 
variables, V932 Sco has otherwise been retained in the sample. 

Table[T]below shows the available information for each of the 
84 stars (V932 Sco, a known Orion variable, was excluded so as to 
provide a list of known and probable LPV stars) and lists the name 
of the known or candidate variable, the name of the star observed 
by STEREO/Hl-l where this is different from the star believed to be 
the origin of the variability, the right ascension and declination (in 
degrees) of the star observed by STEREOfHI-l ,the period found 
from the analysis of the STEREO/HI- 1 lightcurve (in days), the er- 
ror given for this period by PERANSO (in days), which algorithm 
produced this period (1 for PDM, 2 for DFT or 3 for SLM), the 
period given by SlMBAD (NA is shown if none is available) - pref- 
erentially the GCVS ( |Samus et al.|2012| >, whether the GCVS notes 
that the period has been observed to vary, the spectral type given 
by SlMBAD, the variability type given in the GCVS (including the 
New Suspected Variable supplement, these are marked as such), 
where known, followed by the Modified Julian Date (MJD) of all 



maxima observed in the STEREOfHl-1 photometry, a visual esti- 
mate of the error in this time (in days) and the STEREO/HI-1 mag- 
nitude of each maxima. If a spectral type determination has not 
been made or the star has not previously been recorded as variable, 
then " — " is used to denote the absence of this data. 



2 Orion variables are very young stars that have not yet joined the main 
sequence and show large amplitude irregular variability 



Table 1: 84 stars showing long period variability in STEREOMI-l data. V932 Sco was excluded from the table as this 
is a known Orion variable and the new and unclassified variables are more likely to be Miras or semi-regular variables, 
making this table a resource of known and probable LPV stars. 
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Table 2: The dates shown (in MJD-2400000) in this table are when the 6 new candidate variables will next be observed 
by STEREO/HI- 1 . The times take into account the mask that is routinely applied to exclude data likely to be affected by 
solar activity. As the Sun is going to be at maximum during these observations, it is not expected that observations will 
be reliable within the region of the CCDs excluded by the mask. 



star name 


star observed 


RA 


DEC 


start and end dates of upcoming STEREO/HI- ] observations (MJD-2400000) 




by STEREO/HI-\ 


(deg) 


(deg) 


S7EREO/HI-1A 


STEREOMl- IB 


STEREOIHl-lA 


STEREO/m-lB 


SJiTJfiO/HI-lA 


IRAS 16482-2039 


NOMAD 1 0692-0384552 


252.8140 


-020.7219 


56119.30339 


56131.47135 


56418.86914 


56432.20117 


56463.46094-56475.62891 


56807.82813 -56821.16016 


56807.61849 


56819.78646 


IRAS 17289-1917 


NOMAD 1 0706-0429226 


262.9790 


-019.3255 


56129.12500 


56142.40234 


56428.35156 


56441.51953 


56473.65234 - 56486.92969 


56817.29688-56830.46484 


56818.17969 


56831.45703 


IRC -30357 


NOMAD 1 0635-0722217 


271.8500 


-026.4031 


56137.03385 


56150.17057 


56436.99805 


56450.05664 


56481.57031 -56494.70703 


56826.03906 - 56839.09766 


56826.10677 


56839.24349 


IRC -20507 


NOMADI 0714-0715901 


280.0740 


-018.5605 


56144.72005 


56158.02474 


56444.58789 


56457.72852 


56489.23828 - 56502.54297 


56833.54688 - 56846.68750 


56833.75651 


56847.06120 


2MASS J19291709-2034504 


NOMAD 1 0693-0875859 


292.3310 


-020.6250 


56155.52474 


56168.82943 


56456.15234 


56469.20703 


56500.04297 - 56513.34766 


56845.20703 - 56858.26172 


56844.56120 


56857.86589 


NOMAD 1 0784-0674630 




298.3360 


-011.5768 


56162.64714 


56182.14714 


56456.62109 


56476.50781 


56507.15625 -56526.65625 


56845.66406 - 56865.55078 


56851.66536 


56871.16536 
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DISTRIBUTION OF PROPORTIONAL ERRORS (PDM) 



Figure 4. ESO R band image centred on the Mira variable Y Sco (left), 
STEREO/Hl-l A image showing the same field of view near minimum (mid- 
dle) and STEREO/HI- IB image of the same field of view at maximum 
(right). The nearby bright star is NOMAD 1 0706-0340442 and is the only 
star in this field of view in the STEREO database. The variability of Y 
Sco was detected and its period determined by virtue of the blending ef- 
fect with this star, about 2.4 pixels away. The overlaid dark circles illustrate 
the photometric aperture of the Heliospheric Imagers, which is slightly dif- 
ferent for each imager, at 3.2 pixels for STEREO/Hl-IA and 3.1 pixels for 
STEREO/Hl-IB jBewsher et al.|2010) . The galaxy with a bright core just 
below Y Sco has not been recorded. 



The sample of 85 stars for which a period was determined 
were analysed fo r trends in their periods. Some of this analysis was 
done using B^] r Develop ment Core Team] 2008) and some using 
TOPCA!^] The main reason for doing this was that the errors re- 
turned for the periods by PERANSO were often very large and it 
was necessary to verify whether they were valid. Secondly, there 
were no stars in the sample for which the best period found was be- 
tween 200 and 300 days, which was not expected. All these 85 stars 
were also checked for known periods in SlMBAD, the General Cat- 
alogue of Variable Stars (GCVS) (Sa mus et al.|20 12) period being 
used preferentially where multiple determinations were available. 
The actual STEREO/Hl-l images of several stars were examined 
showing maximum and minimum brightness in order to ascertain 
that even from a distance of 2 or 3 pixels, some of the stars in the 
sample could still have been observed indirectly (e.g. Y Sco, shown 
in Figure |4j. 



3.1 Trends in the data 

The concern regarding the accuracy of the errors in the periods re- 
turned by PERANSO is demonstrated by Figures [5] [6] and [7] The 
large size of the errors implies that the algorithms were, individu- 
ally, struggling to resolve a signal, however by comparing the me- 
dian period of the three with the known periods from the GCVS 
l |Samus et al.|2012) shows that a more appropriate la error bar for 
the entire sample is 4% (Figure[8]l. A direct plot of the median pe- 
riod found here against the GCVS periods shows a few cases where 
harmonics may have been found by STEREO/Hl-l instead of the 
correct period but confirms the overall good match (Figure[9]l. The 
table giving all the median periods and other information of rele- 
vance for each star in the sample is found in the appendix in Table 
[JJ Lightcurves phase-folded on these periods are also found in the 
appendix. 

There are no periods found between 200 and 300 days in the 
sample, as shown on the right in Figure [TO] This is in contrast to 
what is expected, as shown by the distribution of the known peri- 
ods of Miras from Kharchenko & Kilpio ( 2000 1 displayed on the 
left in Figure [10] This is believed to be the result of the way the 
sample was selected, with a visual examination of large numbers 



El 



h n n 



3 http://www.R-project.org 

4 http://www.starlink.ac.uk/topcat/ 



0.1 0.2 0.3 0.4 

PROPORTIONAL ERROR IN STEREO PERIOD (PDM ONLY) 

Figure 5. Histogram showing the distribution of the proportional errors in 
the periods observed by STEREO/Hl-l using Phase Dispersion Minimisa- 
tion (Stellingwerf|1978) . 

of lightcurves: periods in this range would have a less distinctive 
pattern in the lightcurve and be almost impossible to distinguish 
from systematic effects. 

The distribution of colours of known Miras from Kha rchenko I 
& Kilpio ( 2000 ) shows two main populations (FigurefTJJ left) when 
comparing the Bl — Rl colours from Mon et et al.| ( |2003) against 
the J — K colours from |Skrutskie et al.| ( |2006| . The distribution of 
colours of stars in this sample observed by STEREO/Hl-l is very 
different, however (Figure[TT] right). This may in part reflect that 
it includes some semi-regular variables but might also indicate that 
some of the stars in the sample have unusual features, such as cir- 
cumstellar dust shells, or that they could have been misclassified. 
Checking the B2 — I colours against J — K, from the same sources, 
shows a population of objects (Figure [T2] left) which is a better 
match for the sample observed by STEREO/Hl-l (Figure[l2] right). 
A slight bias for new and unclassified objects to have large B2 — I 
is expected as a result of checking for Mira-like colours (i.e. very 
red objects) when attempting to ascertain the source of variability, 
however it was not expected that the known Miras found would 
share this feature. The throughput of the STEREO/Hl-l imagers at 
about 950 nm is more likely to be responsible, as these objects are 
bright in the / band. 



4 NOTES ON INDIVIDUAL STARS 

This section details STEREO/Hl-l observations of particular stars 
of individual interest. To clarify the observations of those stars 
without a prior period determination, unfolded lightcurves are 
given for all the new candidates (Figure [T3|>, all the previously un- 
classified variables (Figure [14} and the remaining stars for which 
a period has not previously been determined (Figure \15\. A few 
general comments about the observations are also worth noting: 

• The times of maxima observed do not always correspond to 
the periods given in Table [JJ This is because the times of maxima, 
the shape and amplitude of the lightcurve can all potentially change 
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Figure 6. Histogram showing the distribution of the proportional errors in 
the periods observed by STEREO/Hl-i using a Discrete Fourier Transform 
(Deeming|1975) . 



Figure 8. Histogram showing the distribution in the proportional differ- 
ences ((Pstereo — P<xvs) /Pgcvs) between the STEREO/Hl-l periods and 
those found in the literature (preferentially Samus et al.|2012) . 
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COMPARISON OF PERIODS FROM STEREO AND GCVS 
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Figure 7. Histogram showing the distribution of the proportional errors in 
the periods observed by STEREO/Hl-l using String Length Minimisation 
(Renson|1978) . 



Figure 9. Plot comparing the periods observed by STEREO/Hl-l to those 
found in the literature (preferentially Sam us et al.|2012) . 



from cycle to cycle and the algorithms do not always give a good 
fit. 

• The shape of the lightcurve at maximum can sometimes be 
very flat, such that the star appears constant for many days, even 
longer than a single epoch of observations (about 20 days). In these 
cases, no maximum was recorded as it was not possible to deter- 
mine a central time of maximum with an estimate of the error. 
Therefore, only those maxima sharply defined enough to produce 
an error of about 5 days or less were recorded. 

• In cases where the STEREOfHl-1 period found is a harmonic 



of the known period, it is difficult to know which is genuine, es- 
pecially for those with periods near to a year or a fraction thereof. 
It would be be necessary to combine all available observations to 
be sure for these stars but, although this is beyond the scope of 
the present paper, due to the relatively small number of data points 
herein we are not questioning the accuracy of periods given in the 
literature. 

• Classifications for the new long period variable candidates, or 
previously unclassified variables, cannot be accurately determined 
as the magnitude of variability in the non-standard STEREO/HI-1 
bandpass is not straightforward to convert to a more standard fil- 
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Figure 10. Histograms comparing the periods of known Mira variables 
from Kharchenko & Kilpio 12000), shown on the left, with the distribution 
of periods found for all variables in this work (right). The lack of periods 
between 200 and 300 days in this work is believed to be due to a selection 
effect, with periods in this range being more difficult to distinguish from 
systematic effects. 
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Figure 11. Plot of the Bl-R\ colours (Monet et al.|2003) against J-K 
colours (Skrutskie et al. 20061 for all the Miras in |Kha rchenko & Kilpio] 
(2000) , left, and all the stars in this sample observed by STEREOMl-l, 
right. V2055 Sgr, a known Mira, is the star with the largest value of BY — Ri 
in our sample, whilst V7 1 8 Tau, also a known Mira, has the second-largest 
value of B 1 — R 1 and the largest value of / — K. 



ter. In comparison to the brighter, better-known Miras, however, it 
is likely that most are semi-regular variables, if genuine. The red- 
dening for many of the new candidates, some recorded as having 
circumstellar envelopes, is an additional complication. 

• For 19 stars in the sample, the period presented here consti- 
tutes the first determination of a period. The irregular nature of 
these variables and the large errors in the periods require that these 
are taken with caution, individually, although for those stars with 
a known period there is general good agreement with the periods 
found here and thus there should also be similarly good agree- 
ment, overall, with those for which no period has previously been 
reported. 



4.1 New candidate and unclassified variables 

There are 6 candidate long period variables in the sample not listed 
as variables in the GCVS or its NSV supplement iSam us et al.| 
|2012| > and a further 7 are listed as variables but not classified. Ow- 
ing to the difficulties in determining the amplitude of variability, it 



Figure 12. Plot of the B2-I colours (Monet et al.|20 03 1 against J-K 
colours (Skrutskie et al.|2006) for all the Miras in |Kharchenko & Kilpio| 
(2000) , left, and all the stars in this sample observed by STEREO/Hl-l, 
right. The stars with (B2-I) > 10 from the STEREOMl-\ sample are IK 
Tau, IRC -20507, 2MASS J19291709-2034504, OP Sgr and FQ Sgr, in de- 
creasing order. 



is not possible to ascertain whether these are Miras or semi-regular 
variables, however given their apparent magnitudes and the lack of 
prior measurements or detection of their variability, it is perhaps 
more likely that most are semi-regular variables. In particular, the 
sample was selected partly by the colours shown in Zach arias et| 
[aTl pOO?] ) and are very red objects, thus the likely classifications 
are limited to Miras and semi-regular variables. As five of the new 
candidates were detected indirectly through blending with a nearby 
star, their lightcurves were also individually extracted using the co- 
ordinates provided by SlMBAD, in order to confirm the new sources 
are accurate. The lightcurves of all the new candidates are thus 
presented in Figure [T3] along with a constant star for comparison 
(HD 1651), demonstrating that large amplitude variability has gen- 
uinely been observed in these stars. Note that the constant star has a 
small difference in the observed magnitudes between STEREO/H1- 
1A and STEREO/Hl-IB but that this is on the order of 50 mmag 
and that it remains constant for each satellite. 

The new candidate long period variables are: 

IRAS 16482-2039 There are probable systematic effects confus- 
ing the behaviour near minima. Nevertheless, the amplitude of vari- 
ability clearly exceeds one magnitude. Only the region near mini- 
mum brightness and the decline from maximum are observed, as to 
be expected for a star with a period close to the orbital period of 
the STEREO-Behind satellite. This star is recorded as an infra-red 
source by IRAS (Neugebauer et al. 1984) and has no mention of 
variability. It also appears in the AKAR1 FIS All-sky survey point 
source catalogue I Yamam ura et al.|2010) and has fluxes at 65, 90 
and 140 microns, the flux at 140 microns being equivalent to a mag- 
nitude of about 7.6, with the other two bands registering a magni- 
tude of about 9.9. It was observed to show SiO masing in Deguchi 
|et al.| (2004), which implies that it is an evolved O-rich star under- 
going mass-loss. 

IRAS 17289-1917 This star has either sharp maxima or maxima 
of varying brightness between epochs, although the minima appear 
smooth. The period is very close to one year and any variability 
could thus easily have been difficult or impossible to observe prior 
to the STEREOfHl-1 observations. It is recorded as an IR source by 
IRAS and features inlSkrutskie et al. ( 2006 1 but it does not appear 
in |Yamamura et al7|j2010^ and has not been the subject of further 
investigation. 

IRC -30357 This star shows extremely large amplitude variabil- 
ity but is hampered by the fact that it encroaches on the faint limit 
observable by STEREOfHl-1 . It is possible that only the maxima 
are seen and that the star is too faint to observe at other times. The 
brightness of maximum may be variable between cycles. The star 
has been observed by IRAS and has been assigned a spectral type of 
M8, although there is no mention of variability and it does not ap- 
pear in Yamamura et al. (2010). The STEREOfHl-1 data also shows 
other variability with a period visually near 1 day, probably a result 
of blending with a nearby star possessing rotational variability, al- 
though none are recorded as such. 

IRC -20507 The phase-folded lightcurve is relatively smooth for 
this star but it is unclear if there is a change in the brightness of 
maxima and minima between cycles. One maximum is observed 
but it is possible a second was just caught in the previous cycle, in 
which case the period given here of 43 1 days may be very wrong, 
with the most recent cycle having about 352 days between adjacent 
maxima. This would then be a sign of irregular variability or a sig- 
nificant period change. OH masing has been observed in this star 
l |te LinteiTle kkert et al.|1991[ > and we therefore expect it to have 
a circumstellar envelope, making it a particularly interesting target 
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for future study. The large value of B2 — I ( Monet et al. 2003}, sec- 
ond only to IK Tau, is suggestive of the obscuring presence of a 
circumstellar envelope. The star is designated spectral type M7 in 
SlMBAD. In |Yamamura et al.| < [2"010| > this star has fluxes at 65 and 90 
microns, with equivalent magnitudes of 6.4 and 7.15, respectively. 
Variability with a periodicity visually near 1 day is also seen in 
the data, likely due to blending with a rotational variable, although 
none are recorded as such nearby. 

2MASS J 1929 1709-2034504 No maxima are observed for this 
star but there are indications that the brightness of maximum could 
be variable. The star is listed in SlMBAD as spectral type M7 but 
very few observations have been made and no indication of vari- 
ability is mentioned. The extremely large value of B2 — I ( Monet 
|et al.|20 03 1 is the third-largest in the sample, after IK Tau and IRC 
-20507, which suggests that this star may host a circumstellar en- 
velope. Follow-up observations to check for maser activity would 
be required to confirm this. 

NOMAD 1 0784-0674630 This star shows an erratic unfolded 
lightcurve and the median period of 381 days indicates the shape 
and brightness of maxima and minima may be variable. There were 
no candidate sources near the observed co-ordinates other than this 
star that are likely to be Miras, with two other sources very nearby 
in |Zacharias et al.| ( |2004| l, of which one has colours resembling an 
A-type star and the other a large proper motion typical of a red 
dwarf. Another red star almost 3 arc -minutes away is less likely 
to be the source, TYC 5733-2876-1 (NOMAD 1 0784-0674475), 
owing to its distance, magnitude and the magnitude of variability 
observed, although it should not be excluded from consideration. 

The previously unclassified variables (their unfolded 
lightcurves are shown in Figure [T4|l are: 

CI Sco The STEREO/H1-IA data appears to show a systematic 
trend that may be potentially confusing the period determination, 
although all the STEREO/HI-IB data appear to be at or very near 
minimum, consistent with the period being close to the orbital pe- 
riod of the STEREO-Behind satellite. There is some other variabil- 
ity observed in the lightcurve with a period visually near 5 days, 
although it is not clear in which star this originates. Although the 
variability is unclassified and no period is given, an amplitude of 
1.5 magnitudes has been recorded ( |Samus et al.||2012| >, which is 
about twice that observed by STEREO/HI-1. 

IRAS 16469-3211 Poor phase coverage for this star contributes 
to an uncertainty in both the period and, given the uncertainties in 
the differences in the magnitudes seen by the two satellites, to the 
overall amplitude. One epoch of data clearly shows a sharp increase 
in brightness and the phase-folded lightcurve places this close to 
the maximum in phase. This star is recorded as having an enve- 
lope of OH/IR type as a result of observations reported in |te Lintel] 
|Hekkert et aT] (T99T] > in which OH masing was detected. It was 
previously observed by IRAS but there is no mention of variability. 
The actual source observed by STEREO/HI-1 is nearby (NOMAD 1 
0577-0577145) and has very red colours also and there is a chance 
that this star may be the variable, if it is itself a Mira (however it 
has no known variability of its own or observations other than pho- 
tometry). 

EG Oph For this star the phase-folded lightcurve is smooth ex- 
cept for the first epoch of STEREO/HI-1 A data, which might indi- 
cate a brighter maximum for that cycle. No period or classification 
has been determined for this star but the amplitude is reported as 
1.4 magnitudes (Sa mus et al.|2012 i. It has been observed to be a 
SiO maser source (Deguchi et al. 2004). There is a risk of contami- 
nation from EI Oph, approximately 3 pixels away in STEREO/HI-1, 



which has the same amplitude but is fainter and even more poorly- 
observed. 

HR Sgr No period or classification is given in the GCVS or NSV 
( |Samus et al.|2012| l, however the given amplitude of 4.2 magnitudes 
is suggestive of a Mira. The amplitude observed by STEREO/HI- 1 
of about 1.5 magnitudes is not unusual in the context of the known 
Miras in this sample. There are indications that the signal seen is 
blended with a variable resembling a WUMa type eclipsing binary, 
although this would not be expected to produce variability on the 
scale of a Mira. The median period found here is close to the orbital 
period of the STEREO-Behind satellite but the smoother parts of the 
phase-folded lightcurve are from the STEREO- Ahead satellite so a 
systematic is less likely than changing brightness of the maxima 
between cycles. 

NSV 11552 This star has the lowest amplitude of variability in 
the sample and also one of the shortest periods. The amplitude is 
given in SlMBAD as 1.7 magnitudes, however, so this is unlikely to 
be rotational variability. It is unclear from the STEREO/HI- 1 data 
whether the period might instead be twice that given here. Some 
other variability due to blending is evident in the lightcurve, al- 
though the source is uncertain. The two brightest stars in R nearby 
have significant proper motions indicative of red dwarf stars rather 
than Miras. 

IRAS 19263-1922 The STEREO/HI-1 data is ambiguous regard- 
ing the maximum brightness, which might also be influencing the 
period determination. This star was observed to be a SiO maser 
source in Deguchi et al. ( 2007 1 and was noted for having a highly 
unusual radial velocity relative to the Local Standard of Rest, pos- 
sibly a dynamical effect due to the influence of the Galactic bulge 
bar {Deguchi, Shimoikura & Koike|2010) . This star was observed 
in all four wavelengths by AKARI/F1S, with equivalent magnitudes 
in 65, 90, 140 and 160 microns of about 8.2, 9.7, 8.2 and 7.7, re- 
spectively. Some other variability is evident in the lightcurve, pos- 
sibly due to the nearby variable NSV 12065, this having a period 
visually between 1 and 2 days. 

IRAS 20060-2425 There are indications of changing brightness 
of maxima between cycles for this star. The star has been observed 
by IRAS but is otherwise very poorly observed, although there is 
a possibility that the nearby AN 958.1936, recorded as variable 
in |Luyten| ( |1937| >, could have been an observation of it, or that 
the STEREO/HI-1 observations have detected the variability of this 
star. 
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Figure 13. Unfolded lightcurves of the new candidate long period variables. The data shown here was extracted for their co-ordinates as given in SlMBAD for 
those detected indirectly through the pipeline described in Section |2~3] ( all except NOMAD 1 0784-0674630 ) and a constant star for comparison ( HD 1651 ). 
The clear variability seen here further confirms that such indirect detections are able to recover the large amplitude signals of Miras and semi-regular variables. 
Note that the constant star shows a small discrepancy between the STEREO/Hl-IA and STEREOfHl-lH data but this is not on the scale of variability seen in 
the stars in our sample of long period variables and it also remains constant in each satellite's data, unlike the stars in our sample. 



K. T. Wraight, D. Bewsher, Glenn J. White, W. Nowotny, A. J. Norton and C. Paladini 




TIME ( MJD-2400000 ) 

IRAS 16469-321 1 



TIME ( MJD-2400000 ) 

EG OPH 



TIME ( MJD-2400000 ) 

HR SGR 



10.0 - 
. 10.5 - 



TIME ( MJD-2400000 ) 

NSV - 1552 



TIME ( MJD-2400000 ) 

IRAS 19263-1922 




TIME ( MJD-2400000 ) 

IRAS 20060-2425 




TIME ( MJD-2400000 ) 



Figure 14. Unfolded lightcurves of the unclassified long period variables. 
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Figure 15. Unfolded lightcurves of the previously known variables for which the period found here is the first to be determined. 
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Figure 16. A candidate flare-like event (begins at MJD 2454480) seen in Y 
Sco by STEREO/Hl-IB. 

4.2 Y Sco: flare-like event or starspot? 

A recent search of data from the CoRoT mission found at best a 
tentative sign of a single flare-like event |Lebzelter 20111 in only 
one star. The larger sample we present here has the difficulty that 
many potential flares could easily be mistaken for a de-pointing 
event caused by a micrometeorite impact. Y Sco, however, displays 
an event that does not look like any of the known systematics and 
therefore might be a flare, although the timescale of this feature, 
being of about 1 day duration, is longer than might be expected 
(Figure [16}. The feature begins in the first epoch of STEREOIHl- 
1B data, at MJD 2454480 and lasts for one day, during which the 
brightness increases in a linear fashion by 0.1 magnitudes. There 
is no change in the scatter of the lightcurve during this time, thus 
it is not due to a micrometeorite hit. An alternative explanation for 
this feature might be that it is due to a dark spot on the surface 
rotating out of view, or perhaps a bright spot rotating into view. 
The star is only a few days away from maximum brightness at the 
time this event occurs. There is a remark in the GCVS ( |Samus et| 
|al.|2012[ ) that this star has been irregular since 1972 and although 
the period observed by STEREO/Hl-l of 355 days is very similar 
to the pre- 1972 period of 351.88 days, the interval between the two 
observed maxima is about 782 days. If the flare-like event is really 
due to a starspot, it might also mean that the maximum observed a 
few days before this was not a full maximum, and the star was still 
approaching maximum brightness. 

Some caution needs to be taken with the variability, as the star 
actually observed by STEREO/Hl-l is NOMAD 1 0706-0340442, 
2.4 pixels away. There is very vague short period variability with a 
period visually close to a day in the lightcurve, with an amplitude of 
<50 mmag that may be from this star and its colours in |Zacharias| 
|et al.|p004} are suggestive of a star of spectral type G or K. It thus 
has some potential for rotational variability and flares and spots of 
its own, although no other indication of such an event is seen in the 
lightcurve. If it has a rotational period near 1 day, this would likely 
exclude it as the source of the flare-like event, since this lasts about 
as long as one such rotation and could not be due to an event on a 
part of the surface of this star rotating into or out of view. There is 
an unrecorded galaxy with a bright core very close to Y Sco (RA: 
247.356 degrees DEC: -19.3638 degrees) but this is probably too 
faint for any AGN-type effects to be observed by STEREO/Hl-l. 

4.3 Period Changing Variables 

15 of the stars in the sample are recorded in Samus et al. (2012) as 
showing variable periods. Column 9 of Table[T]labels these with a 
1. These stars are briefly discussed here for indications of further 
period changes or other unusual features. Below this list two other 
stars showing evidence of period changes are detailed. 

V Tau This star has the shortest period of any found in the sample, 



at 170 days. This is a very good match to the known period of 
168.7 days, thus no sign of period changing is observed. 

Z Tau This star has the second-longest previously known period 
of any in the sample, at 466.2 days, whilst the period found here is 
just outside the typical 4 % difference, at 438 days. The accuracy of 
the period determination is unable to exclude a small period change 
but a significant change is unlikely. 

U Ori The period found here of 372 days is a good match to the 
known period of 368.3 days. The period is probably constant for 
the duration of the STEREO/Hl-l observations. 

T CMi The period found here of 325 days is a very good match 
to the known period of 328.3 days. A change in the shape between 
cycles is suggested by the phase-folded lightcurve, although a pe- 
riod change might also manifest in this way and not be picked up, 
thus it is inconclusive for this star whether or not any change has 
occurred. 

S Leo The period found here of 197 days is in fair agreement 
with the known period of 190.16 days. There is evidence of a small 
change in the brightness of maxima between cycles. Two succes- 
sive maxima are observed with an interval of 186.3 days between 
them, within the errors of the times of maxima no period change is 
therefore found. 

SS Vir The period found here of 377 days is within the typical 4 % 
away from the known period of 364.14 days. Two successive max- 
ima are observed with 342. 14 days between them, suggesting a pe- 
riod decrease, however. The phase-folded lightcurve is not smooth 
and a period change or irregularity could have resulted in a poor 
period determination. 

S Vir The STEREO/Hl-l period of 356 days is slightly different 
to the 375.1 days from Samus et al. (2012)andthe poor phase cov- 
erage may be hampering the period determination. It is therefore 
uncertain if the period is changing, however the brightness of both 
maxima and minima might also be changing, which would further 
complicate the analysis. 

S Lib The period from STEREO/HI- 1 of 197 days is in good 
agreement with the known period of 192.9 days and it seems un- 
likely that significant period changes have occurred. The phase- 
folded lightcurve nevertheless clearly shows that the brightness of 
maxima and possibly the shape of the lightcurve may differ from 
cycle to cycle. 

Y Sco See above discussion in Section l4~2l 

R Oph In spite of the large uncertainty in the period, the 
STEREO/Hl-l period of 302 days is in good agreement with the 
known period of 306.5 days. There is an indication of the brightness 
of minima varying between cycles but no indication of a significant 
period change. 

VI 869 Sgr The STEREO/Hl-l period of 318 days only moder- 
ately agrees with the known period of 332 days and the phase cov- 
erage is poor as there is less data than normal for this star, likely 
due to it being so near to the edge of the field of view it may have 
been missed completely on some orbits. The analysis is therefore 
unable to determine if there have been any period changes. 

V3876 Sgr The period found here of 344 days is in good agree- 
ment with the known period of 352 days. A period change therefore 
seems unlikely. There is some indication of a change in maximum 
brightness between cycles. 

RX Sgr The period found here of 326 days is in reasonable agree- 
ment with the known period of 335.23 days. An additional compli- 
cation here is the presence of the known semi-regular variable BH 
Sgr about 95 arc-seconds away. Changes in maximum brightness 
between cycles indicated by the lightcurve are therefore unreliable 
and it is uncertain what effect this other variable would have on the 



observations of the period, although with a period from |Samus et| 
[aT|j2012) of 100 days and being a much fainter object, the effects 
of RX Sgr should dominate. 

AN Sgr The period found here of 325 days agrees moderately 
well with the known period of 337.56 days. The phase coverage is 
mostly concentrated around minimum and there is therefore no ev- 
idence of period changing or other significant differences between 
cycles. 

RR Sgr The period found here of 335 days is in excellent agree- 
ment with the known period of 336.33 days. The period is there- 
fore unlikely to be significantly changing, however the unfolded 
lightcurve indicates differences in maximum brightness and possi- 
bly shape between cycles. 

Two other stars show evidence of a different period from that 
given in |Samus et al.| ( |2012) . Both are classified as semi-regular 
variables rather than Miras. Their unfolded lightcurves are shown 
in Figure[l7] Unfortunately, both are too poorly-observed from the 
ground for the observations to be confirmed - no data for them 
exists in the online archives of the AAVSO, the BAA VSS or 
the AFOEV. Their position on the Ecliptic Plane makes long-term 
monitoring difficult and without long-term monitoring it is impos- 
sible to ascertain whether the periods are actively changing or me- 
andering about a more regular value (Zijlstra & Bedding 2002 1. The 
two stars are 

V5545 Sgr Although the period found here of 368 days is in rea- 
sonable agreement with the known period of 377 days, the three 
consecutive times of maximum brightness observed indicate a pe- 
riod of about 388 days (Figure [T7] upper plot). A small period 
change may therefore have occurred. 

V360 Sgr The period found here of 367 days is, at best, a 
very rough approximation of a harmonic of the known period of 
165 days. The three maxima observed are not consistent with ei- 
ther period, being separated by 383.47 days and 389.39 days (Fig- 
ure|17| lower plot). A significant period change may therefore have 
occurred, or the period may be irregular. 



5 DISCUSSION 

The arrangement of two almost identical cameras on almost identi- 
cal satellites in different heliocentric orbits permits the observation 
of periodicity on timescales not observable from the Earth or from 
a single satellite alone. The unusual bandpass of the STEREOfHI-1 
instruments, in particular the throughput in the infra-red at about 
950 nm (Figure [TJ, allows for the observation of stars of very late 
spectral type and also for stars that are obscured by circumstel- 
lar dust shells. Together, this provides homogenous observations of 
Mira variables with improved phase coverage and has potentially 
detected some new objects of this type or the related semi-regular 
class. 

Even though some objects have either no R magnitude listed 
in |Zach arias et al. (2004) or are fainter than the cutoff of 12th mag- 
nitude used to select objects to be observed by STEREOfHI-l, the 
large pixel scale of 70 arc-seconds nevertheless allows them to have 
an indirect influence on the lightcurves of very nearby stars that 
are in the STEREO/Hl-1 database (Figure 01. Although in some 
regions, the Galactic centre especially, it is difficult or impossible 
to ascertain the likely source of variability, the very large ampli- 
tudes and very long periods of Miras and semi-regular variables is 
sufficiently distinctive that in some cases the variability is reliably 
observed even in a densely populated field of view. In the case of 
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Figure 17. The unfolded STEREOMl- 1 lightcurves of V5545 Sgr (upper) 
and V360 Sgr (lower). Maxima for V5545 Sgr are observed at (times are 
MJD-2400000): 54506.99 ±2.5; 54895.71 ±3.0 and 55283.38 ±3.0. Max- 
ima for V360 Sgr are observed at (times are MJD-2400000): 54506.99 
±2.5; 54895.71 ±3.0 and 55283.38 ±3.0. 



new candidate variables, it further complicates the matter of clas- 
sification, which is essentially impossible from STEREO/Hl-l data 
alone, although it is likely most are of the semi-regular class. 

The majority of the sample of 85 stars are known Mira vari- 
ables, with a small minority classified as semi-regular variables and 
a single variable of the Orion type. 7 are listed as variable in either 
the GCVS or NSV |Samus et al] ( |2012) but without a classifica- 
tion and 6 have not been previously observed to be variable. For 
19 stars, the period we present here is the first determination of a 
period. There is a reasonable agreement in the periods previously 
known and those found by an analysis of the STEREOfHl-l pho- 
tometry (Figure|9]l, although a harmonic is occasionally found and 
also some of the stars are known to change their periods. The ac- 
curacy of the period determination is insufficient to detect small 
changes in period of a few days but larger changes of a couple 
of weeks are potentially observable (Figure [8j. The new candidate 
and unclassified variables are individually discussed, as are those 
known or suspected of period changing. One star showing a partic- 
ularly unusual feature is discussed: Y Sco, which shows a candidate 
flare-like event, although a starspot might also be an alternative ex- 
planation (Figure [T6j. 

Of the newly-discovered variables, two (IRC -20507 and 
2MASS J19291709-2034504) have exceptionally large values of 
B2-7 (Figure \12\. This may be suggestive of the presence of a 
circumstellar envelope and these stars might be therefore poten- 
tial maser sources, indeed IRC -20507 has already been recorded 
as such (te Lintel Hek kert et al.| 19 91 ). Follow-up observations of 
these stars would contribute to the understanding of the evolution 
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of stars during the AGB stage. There are features of some stars of 
this type that cannot be explained by current models, such as an 
excess of ammonia (Menten et al. 20101 and sub-solar values of 
16 0/ n O and 16 Q/ 18 Q (jDecin et al.feoTof, and having more exam- 



ples to study would be very useful. We encourage interested readers 
to make follow-up observations of the new candidate variables, in 
order to confirm their nature and periodicity. To assist in this, we 
provide dates when the new candidate variables will be observed 
by STEREOmi-i so that observations may be conducted simulta- 
neously (Table[2](. It is hoped that five years' worth of STEREO/H1- 
1 data will soon be made available for all stars in the field of view 
with listed R magnitudes of 12 or greater in the NOMAD 1 cat- 
alogue (Zacharias et al. 2004), although a date has not yet been 
fixed. In due course, more data will be gathered by STEREO/HI but 
it is not known when, or even if, this will be made available on a 
large scale. 

A complete lack of objects with periods between 200 and 
300 days is observed in the sample (Figure [To] right). There are 
a small number found with periods that are over 400 days thus it 
cannot be completely excluded that one or more of those may in 
fact be harmonics with a genuine period in this range. Nevertheless, 
for a sample of this size, this is a significant feature and the only 
likely explanation is that it is the result of a selection effect. In the 
earliest stage of the analysis, many thousands of lightcurves were 
visually examined in a search for long period variability, however it 
is not easy to distinguish artificial effects from genuine variability, 
especially for stars with potentially very erratic variability in terms 
of the magnitude of maxima, minima, shape and even period. This 
may have caused genuine variables with periodicity in this range to 
be too unconvincing to have been recorded as a likely variable. A 
similar selection effect might also have led to a preference for find- 
ing variables with periods near 1 year or close to the orbital periods 
of the two STEREO satellites as these produce a more recognisable 
pattern in the unfolded lightcurves. It is therefore likely that more 
Miras and semi-regular variables remain undiscovered in the data. 

A selection effect may also be responsible for the prevalence 
of stars in our sample with large B2 — I values from |Monet et] 
[aT] |2003|, shown on the right in Figure [12] The search for the 
sources of variability was focused on very red objects, although 
other colours are also given in Zacharias et al. (2004), in particular 
the J, H and K colours from Skrutsk ie et al.| p006l. This preva- 
lence could also be partly due to the sensitivity of the STEREO/H1- 
1 imagers at about 950 nm (Figure [TJ, making it easier to recover 
variability from stars bright at this wavelength. 

The photometry gathered by STEREO/HI- 1 on both long and 
short timescales of bright stars on the Ecliptic Plane is a valuable 
resource for the observation and monitoring of Miras and semi- 
regular variables. Some variables of this type have periods close to 
a year, or a fraction thereof, and the full range of their behaviour 
cannot be monitored from Earth or Earth orbit and thus the ho- 
mogenous observations of STEREO/Hl-l provide a unique window 
to advance the study of these poorly-understood objects. 
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APPENDIX A: PHASE-FOLDED LIGHTCURVES 

The phase-folded lightcurves for all 85 stars in the sample are 
shown in Figures [AT| | A2| |A3] | A4| and | A5 1 In each case, the period 
folded on is shown in the title, along with some other information 
on the star for ease of reference. 



ACKNOWLEDGEMENTS 

The Heliospheric Imager (HI) instrument was developed by a 
collaboration that included the Rutherford Appleton Laboratory 
and the University of Birmingham, both in the United Kingdom, 
and the Centre Spatial de Liege (CSL), Belgium, and the US 
Naval Research Laboratory (NRL), Washington DC, USA. The 
STEREO/SECCHI project is an international consortium of the 



STEREO observations of long period variables 



RT (Jam P=J4E. LW5 C Him 



LI Gem P=3B5 LW1S VBc SR 



I Gem P=5B7 DAIS MID Mira 




PH*3E 



PH*SE 

PCnc P=179 OivS Mollis Mira 



w Cnc P=M9 Dai? M7s Mira 




4,4 

s " 




\ ^\ ; 




: \ \ 


I 50 






a 7. Q 


- \ \\ 


I sa 




\ \ - 


I 


> \ \ 


1 ^ 




r"' . r 


| EkO 


y/ . V/ . 


□J 0.5 


1.D 1.E- 2 
PH-5E 


□ Q£ OS IjQ 1.5 2 
FH"5E 



Z Cut P-16S DavS H5« Mfra 



T vlr P-330 [inS H6b Hlra 







\ 









^ - 


s 


I 1 


I } 








l 






£ 9,5 



PWSE 

SS * P-37-7 LWS C SB 



PHASE 

S Pss P-413 LWVS M7e Hlra 



PH»SE 

S ulr P-3&S DAYS M7llla Mini 



1 3.4 

B i.B 

T =■» 

I 

s B - 2 
is 

t5 H-* 





■ ■ ^ 




TJO 




v\ 


5 






- V 




-i 


*> ■ i 










| tui 




/ 


/ i 












t 


f \ 


















. 1 








s 


1 





i 


r 




■_ 


7 


''! 








i 













0.5 1J) 15 

PH*5E 

K.5 Lib P-362 LW5 Me Mini 



EP Lib F— 131 DATT 



1.5 
Mlra 



UD 05 Ml 1.5 

PHASE 

5X Lib F-524 DATE W5c Mlna 






\ 












u 



















V 

\ ': 








Figure Al. Phase-folded lightcurves of stars showing long period variability in STEREO/HI- 1 . 2 phases of the periods indicated are shown. 
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Figure A2. Phase-folded lightcurves of stars showing long period variability in STEREO/HI- 1 . 2 phases of the periods indicated are shown. 
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Figure A3. Phase-folded lightcurves of stars showing long period variability in STEREOfRi- 1 . 2 phases of the periods indicated are shown. 
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Figure A4. Phase-folded lightcurves of stars showing long period variability in STEREO/HI- 1 . 2 phases of the periods indicated are shown. 
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Figure A5. Phase-folded lightcurves of stars showing long period variability in STEREOfRi- 1 . 2 phases of the periods indicated are shown. 
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